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The theoretical status of squark and gluino production at present and future hadron 
colliders is reviewed. I shall concentrate on the evaluation of SUSY-QCD correc- 
tions to the production cross sections and their phenomcnological implications. 

1 Introduction 

Supersymmetry is a new symmetry between fermions and bosons, which has 
not been observed so far at the level of elementary particles. The symmetry 
is motivated by the interest of providing a natural solution to the hierarchy 
problem due to the absence of quadratic divergences in perturbation theory, 
in contrast with the Standard Model. Supersymmetric grand unified theo- 
ries have predicted a Weinberg angle in striking agreement with the present 
measurements at the electroweak scale. An important constraint of super- 
symmetric theories is that they predict a light scalar Higgs boson with mass 
M h < 130 GeV. 

The search for Higgs bosons and supersymmetric particles provides a 
strong motivation for present and future experiments. The colored super- 
symmetric particles, squarks and gluinos, can be searched for at the present 
Tevatron experiments, a pp collider with a cm. energy of 1.8 TeV upgraded 
soon to 2 TeV, and the LHC in the near future, a pp collider with a cm. energy 
of 14 TeV. The Tevatron searches set presently the most stringent mass bounds 
on these particles: at 95% CL, gluinos have to be heavier than ~ 175 GeV, 
while squarks with masses below ~ 175 GeV have been excluded for gluino 
masses rrig < 300 GeVli. In i?-parity conserving theories, e.g. the minimal 
supersymmetric extension of the Standard Model [MSSM], supersymmetric 
particles can only be produced in pairs. All supersymmetric particles will de- 
cay to the lightest one [the LSP], which is stable stable thanks to conserved 
i?-parity. As a consequence the typical signatures for the production of super- 
symmetric particles will mainly be jets and missing transverse energy, carried 
away by the LSP. 

Squarks and gluinos can be produced in 4 different processes, 

pp/pp^ qq,qq,gg,qg + X (1) 

We worked in the approximation of mass-degenerate light-flavored squarks, so 
that mixing effects are absent. However, this approximation is not valid for 
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stop production, where mass splitting and mixing effects are sizeable. Stop 
production has been presented elsewhere □. The light quarks (u, d, s, c,b) will 
be treated as massless particles, while the top quark will be included in loops 
with a mass mt = 175 GeV. Moreover, stops are taken into account in the 
virtual loops, too, with a mass equal to the other squark masses. Since their 
effect will be suppressed, this provides a reasonable approximation. In Eq. [j] 
a summation over all possible squark flavors and charge conjugate final states 
should be implicitly understood. 



2 Lowest Order 



The calculation of the lowest order n (LO) cross sections of the processes in 
Eq. [j] was performed a long time agoS. The 4 production processes divide into 
several partonic subprocesses [Np = 5]: 
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mi = m? — m~ A = y (s — m? — m?) 2 — 4m?m? 

The LO hadronic cross sections can be obtained by convoluting the partonic 
ones with the corresponding parton densities. At the Tevatron qq production is 
the dominant process for squark masses lighter than the gluino mass, while in 
the reverse case gluino pair production dominates. The associated qg produc- 
tion process always provides a sizeable fraction of SUSY particle production. 
Beyond the excluded mass range the qq intial state dominates over the gg ini- 
tial state in qq and gg pair production. At the LHC the situation is reversed: 
for smaller squark masses, qq and qg production are dominant, while for larger 
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squark masses SUSY particle production is dominated by gg and qg produc- 
tion. Moreover, the gg initial state clearly dominates gluino pair production. 
The LO scale dependence leads to an estimate of ~ 50% uncertainty for the 
production cross sections. Since the reconstruction of squarks and gluinos will 
be very difficult due to the escaping LSP, an important possibility to mea- 
sure the squark and gluino masses is provided by the values of the total cross 
sections. Thus a complete NLO calculation is needed in order to reduce the 
theoretical errors to a reliable lcvela. 



3 SUSY QCD Corrections 

The evaluation of the SUSY QCD corrections will be exemplified for qq pro- 
duction. These corrections split into two pieces, the virtual ones, generated by 
virtual particle exchanges, and the real ones, which originate from gluon radi- 
ation and the corresponding crossed processes with three-particle final states. 

3.1 Virtual corrections 



g 'Tnrs~5~5> 
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Figure 1: Typical diagrams of the virtual corrections. 



The one-loop virtual corrections are built up by gluon, gluino, quark and 
squark exchange contributions [see Fig. pj]. Their contraction with the LO 
matrix elements provides the one-loop contributions to the physical matrix 
elements. The evaluation of the virtual corrections has been performed in 
dimensional regularization, leading to the extraction of ultraviolet, infrared and 
collincar singularities as poles in e = (4 — n)/2. For the chiral 75 coupling we 
have used the naive scheme, which is justified in the present analysis at the one- 
loop level without anomalies. We have explicitly checked that after summing 
all virtual corrections no quadratic divergences are left over, in accordance 
with the general property of supersymmetric theories. The renormalization 
has been performed by identifying the squark and gluino masses with their 
pole masses, and defining the strong coupling in the MS scheme, including 
five light flavors in the corresponding function. The massive particles, i.e. 
squarks, gluinos and top quarks, have been .-decoupled by subtracting their 
contribution at vanishing momentum transferH. In dimensional regularization, 
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there is a mismatch between the gluonic degrees of freedom [d.o.f. = n — 2] and 
those of the gluino [d.o.f. = 2], so that SUSY is explicitly broken. In order to 
restore SUSY in the physical observables in the massless limit, an additional 
^jnite counter-term is required for the renormalization of the novel qgq vertex 

3.2 Real corrections 




Figure 2: Typical diagrams of the real corrections. 



The real corrections originate from the radiation of a gluon in all possible 
ways and from the crossed processes by interchanging the gluon of the final 
state against a light quark in the initial state. The phase-space integration of 
the final-state particles has been performed in n = 4 — 2e dimensions, leading 
to the extraction of infrared and collinear singularities as poles in e. In order 
to isolate the singularities we have introduced a cutoff A in the invariant 
mass of, say, the qg pair, which separates soft and hard gluon radiation. After 
evaluating all angular integrals and adding the virtual and real corrections, the 
infrared singularities cancel. The left-over collinear singularities are universal 
and are absorbed in the renormalization of the parton densities at NLO. We 
defined the parton densities in the conventional MS scheme including five light 
flavors, i.e. the squark, gluino and top quark contributions are not included 
in the mass factorization. Finally we end up with an ultraviolet, infrared and 
collinear finite partonic cross section, which is independent of the cutoff for 
A^O. 

However, there is an additional class of physical singularities, which have 
to be regularized!! In the second diagram of Fig. |^ an intermediate qg* state 
is produced, before the [off-shell] gluino splits into a qq pair. If the gluino 
mass is larger than the common squark mass, and the partonic cm. energy is 
larger than the sum of the squark and gluino masses, the intermediate gluino 
can be produced on the mass shell. Thus the real corrections to qq production 
contain a contribution of qg production. The residue of this part corresponds 
to qg production with the subsequent gluino decay g — > qq, which is already 
contained in the LO cross section of qg pair production, including all final-state 
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cascade decays, 



da res (gq -> gq -> qqq) 



dpi 



°(gq -» gq)BR(g ->■ g<?) 



rrigTg/w 



(^_ m 2 )2+m 2 r 2 



(9) 



Thus this term has to be subtracted in order to derive a well-defined production 
cross section for each individual final state. Analogous subtractions emerge 
also in other reactions: if the gluino mass is larger than the squark mass, the 
contributions from g — > qq, qq have to be subtracted, and in the reverse case 
the contributions of squark decays into gluinos have to subtracted. 
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Figure 3: Scaling functions of the gluino pair production cross sections [rj = s/(4m?) — 1/. 
Top mass: rat = 175 GeV. 
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4 Results 



4-1 Partonic cross sections 

The partonic scaling functions are defined as 



a 2 s (Q 2 ) 



fZ +S + /«+/« log" 



12 



(10) 



The parameter m denotes the average mass of the produced massive particles. 
The scaling functions are shown for gluino pair production in Fig. [3] as a func- 
tion of r) = s/ (4m~) — 1 for the gg, qq and gq initial states. The singularity of 
f gq at s = (mq+mg) 2 Jis a left-over singularity of the residue subtraction, which 
however is integrableQ. The scaling functions can be evaluated analytically in 
the threshold region (3 <C 1 and the high energy limit s 3> rn~, m~. 

In the threshold region the dominant contributions originate from Coulomb 
singularities in the f^ +S functions due to Coulombic gluon exchange between 
the final state particles and from large logarithmic enhancements in the flf and 
fij functions due to soft gluon radiation from the initial states. The leading 
logarithms of soft gluon radiation can be resummedQ. The analytic results in 
the threshold region for e.g. gg — > qq are given by 

7 , T ^ »T/ I C »R 11 



192 ' Jga J 99 336/3 

(ii) 

In the high energy region the LO cross sections scale as &lo °c o> s /s, 
while the NLO cross sections develop a high-energy plateau, 3-nlo °c a 2 /m 2 , 
originating from t-channel gluon exchange between an initial quark/gluon line 
and the hard process. Thus the scaling functions fij approach finite values, 
which, can be calculated by using the factorization in kx of the exchanged 
gluoifi. For qq pair production the non-zero limits read as 

H= 2W9_ H= 2159 

J aa jnon_ J 00 
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4-2 Hadronic cross sections 
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Figure 4: K factors of the different squark and gluino production cross sections at the 
Tevatron. Parton density: GRV(94) with Q = m. Top mass: mt = 175 GeV. 
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Figure 5: K factors of the different squark and gluino production cross sections at the LHC. 
Parton density: GRV(94) with Q = m. Top mass: mt = 175 GeV. 



The hadronic cross sections can be obtained by convoluting the partonic 
ones with the corresponding parton densities. We have performed the nu- 
merical analysis for the Tevatron and the LHC. For the natural renormaliza- 
tion/factorization scale choice Q = m, where m denotes the average mass of 
the final-state SUSY particles, the SUSY QCD corrections are large and posi- 
tive, increasing the total cross sections by 10-90%U. This is shown in Figs. H,|l 
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Figure 6: Scale and partem density dependence of the total gg and gq production cross 
sections at the Tevatron in LO and NLO. Parton densities: GRV(94) (solid), CTEQ3 



(dashed) and MRS(A') (dotted); mass parameters: 
m t = 175 GeV. 



280 Ge V, m~ q = 200 Ge V and 



where the K factors, defined as the ratios of the NLO and LO cross sections, 
are presented as a function of the corresponding SUSY particle mass for the 
Tevatron and the LHC. 

We have investigated the residual scale dependence in LO and NLO, which 
is presented in Fig. ||. Including the NLO corrections reduces the LO scale 
dependence by a factor of 3-4 and reaches a typical level of '-s-,15%, which 
serves as an estimate of the remaining theoretical uncertainties a Moreover, 
the dependence on different sets of paxton densities is rather weak and leads 
to an additional uncertainty of ~ 10% In order to quantify the effect of the 
NLO corrections on the search for squarks and gluinos at hadron colliders, we 
have extracted the SUSY particle masses corresponding to several fixed values 
of the production cross sections. These masses are increased by 10-30 GeV 
at the Tevatron [see Fig. and by 10-50 GeV at the LHC [see Fig. §, thus 
enhancing the present bounds on the squark and gluino masses significantlyQ. 



Finally we have evaluated the QCD-corrected transverse-momentum and 
rapidity distributions for all different processes. As can be inferred from Fig. pi 
the modification of the normalized distributions in NLO compared to LO is 
less than 10% for the transverse-momentum distributions and negligible for 
the rapidity distributions. Thus it is a sufficient approximation to rescale the 
LO distributions uniformly by the K factors of the total cross sectionsu. 
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Figure 7: The total cross sections for the Tevatron [y/S =1.8 TeV]. NLO (solid): GRV(94) 
parton densities, with scale Q = m; compared with LO (dashed): EHLQ parton densities, 
at the scale Q = y/s. 



5 Conclusions 



We have reviewed the theoretical status of the squark and gluino production 
cross sections at present and future hadron colliders. The evaluation of the full 
SUSY-QCD corrections has been described in detail. They enhance the cross 
sections by about 10-90% and are thus significant for a reliable prediction. In 
spite of the large size of the NLO corrections, the residual scale dependence 
reduced by a factor of 2.5-4, so that the theoretical uncertainty decreased 
to ~ 15%. The dependence on different parton densities is ~ 10% and thus 
weak. The NLO corrections increase the bounds on the squark and gluino 
masses, extracted from upper bounds on the total cross sections, by about 
15-35 GeV at the Tevatron and by 10-50 GeV at the LHC. The shape of the 
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Figure 8: The total cross sections for the LHC [yS = 14 TeVJ. NLO (solid) compared with 
LO (dashed). GRV(94) parton densities, with scale Q = m. 

differential distributions in the transverse momentum px and the rapidity y are 
hardly affected by the NLO corrections, so that a simple rescaling of the LO 
distributions by the K factors of the total cross sections provides a reasonable 
approximation within^ 10%. The final results are encoded in the Fortran 
program PROSPINol. f] 
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